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ABSTRACT 

Stellar winds are a crucial component of massive stars, but their exact properties still remain uncertain. 
To shed some light on this subject, we have analyzed an exceptional set of X-ray observations of f Puppis, 
one of the closest and brightest massive stars. The sensitive lightcurves that were derived reveal two major 
results. On the one hand, a slow modulation of the X-ray flux (with a relative amplitude of up to 15% over 
16h in the 0.3^.0keV band) is detected. Its characteristic timescale cannot be determined with precision, 
but amounts from one to several days. It could be related to corotating interaction regions, known to exist 
in ^Puppis from UV observations. Hour-long changes, linked to flares or to the pulsation activity, are 
not observed in the last decade covered by the XMM observations; the 17h tentative period, previously 
reported in a ROSAT analysis, is not confirmed either and is thus transient, at best. On the other hand, 
short-term changes are surprisingly small (< 1% relative amplitude for the total energy band). In fact, 
they are compatible solely with the presence of Poisson noise in the data. This surprisingly low level 
of short-term variability, in view of the embedded wind-shock origin, requires a very high fragmentation 
of the stellar wind, for both absorbing and emitting features (> 10^ parcels, comparing with a 2D wind 
model). This is the first time that constraints have been placed on the number of clumps in an O-type star 
wind and from X-ray observations. 

Subject headings: X-rays: stars - Stars: early-type - Stars: individuals: f Puppis 



1. Introduction 

Stars more massive than 20 suns are the true 'cos- 
mic engines' of our Universe. Through their winds and 
final explosions as supernovae, they shape the inter- 
stellar medium and largely contribute to its chemical 
enrichment. Besides, those stars are the most lumi- 
nous ones, the only ones seen from afar. Studying dis- 
tant galaxies or our close neighborhood thus requires 
a good knowledge of these objects. However, one of 
their most crucial physical properties, the stellar wind. 



remains poorly constrained. 

In recent years a large debate arose in the massive 
star's community regarding the structure and strength 
of such winds. Based on the a nalysis of resonance 



lines in UV and FUV spe ctra, iBouret et al 



and FuUerton et al 



(I2005h 

mass-loss rates 



(120061) revised 
of a sample of O-type stars by orders of magnitude 
down compared to those obtained from Ha and ra- 
dio diagnostics. lOskinova et al.l (120071) showed that 
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including wind clumping in the analysis of UV reso 
nance lines yielded mass-loss r ates in agreement with 
those derived from Ho- (see also lSundqvist et al I2OIII 



Surlan et alj|2012l) . These estimates crucially depend 



on the degree of wind clumping and exact properties 
of these clumps. 

In this debate. X-ray observations may well play 
a crucial role. Indeed, the first high-resolution spec- 
tra of massive stars revealed line profiles which were 
more symmetric and le ss blueshifted than expected 
(e.g. iKramer et alj|2003h . triggering (for a large part) 
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the current mass-loss debate. The first global fitting of 
high-resolution X-ray spectra has just been attempted 
and it is showing the clumps to be rather spherical (and 
not "pancake-like"), the non-porosity of the wind, and 
some of the plasma emitting components to be dis- 
tributed over large regions ( up to a hundred stellar 
radii in size, [Herve et al.ll2012[ and Herve et al., sub- 
mitted). 

Another observable parameter, the variability, 
could also be used to learn about the wind structure. 
Indeed, the intrinsic instability of line-driven stellar 
winds generates clumps and shocks between clumps, 
which in turn produce high - energ y emission over 



a large zone (iLucy & Whitd 119801 : iFeldmeier et al. 
1997bh . This "embedded wind-shock model" cur- 
rently is the standard model for X-ray production. Be- 
cause the X-rays are directly born in the winds, their 
variability provides strong constraints on the wind 
structure. Hydrodynamical simulations showed the 
highly changing level of the produced X-rays, a vari- 
ability which could only be lowered if many clumps 



were present (IFeldmeier et al. I ll997bh . A quantitative, 
theoretical assessment of the interplay between clump 
number and variation level, for different energ y rang es. 



number and variation level, to r dilierent energ 
was performed only recently jOskinova et alj 



2004 . 



On the observational side, searches for such high- 
energy variations were undertaken with several fa- 
cilities, the m ost thorough work having be en done 
with ROSAT (iBerghofer & SchmittI 1 1994b). but no 
significant short-term variation was detected in O- 
stars. However, the sensitivity of these early obser- 
vations was poor, with relative error bars in ROSAT 
Ughtcurves amounting to several tens of percents, and 
they thus only provided weak constraints. 

In the last decade, XMM-Newton has obtained an 
exceptional IMs dataset of a massive star, ^Puppis. To 
date, there exists no more sensitive dataset, or a dataset 
with a better (time and spectral) coverage, for a mas- 
sive star. We have thus undertaken the analysis of this 
dataset, with the aim of constraining the wind struc 



ture u sing X-ray variability. A first paper (iNaze et al. 



20ll hereafter Paper I) has presented these data, ex- 
plained their reduction, and provided some first results 
(notably the analysis of EPIC spectra). This second 
paper focuses on the detailed variations recorded with 
EPIC and RGS instruments. It is organized as follows: 
the target is presented in Sect. 2, and the observations 
in Sect. 3; the observed changes found on long, inter- 
mediate, and short timescales are described in Sects. 
4 to 6; the results are finally summarized in Sect. 7. 



An appendix provides details on the results for each 
exposure. 

2. ^Puppis as a variable star 

Being one of the closest and brightest massive stars, 
^ Puppis (04Infp) provides the opportunity to perform 
an in-depth investigation of its properties, notably its 
variability. 

For example, ( Puppis is well known for its spectral 
variations in the visible and UV domains. Some im- 
pressive variations have been report ed, e.g. a 50% de- 
creas e of Ha strength over 2 years dConti & Niemelal 
1976h and a doubhng of the Hen .14686 and Ha 
line profiles over successive nights (i.e. from si n- 



gle li ne to double-peaked profiles, IWegner & Snow 
1978h . Periodic changes have also been detected, 
on short as well as longer timescales. A 5d mod- 
ulation found in the UV and optical data (Ha line. 
Si IV doublet at 1400A, photometry in the Stromgren b 
band) was interpreted as linked to the stellar rotation 



(Moffat & Michaud 1981; Balona 1992; Howart h et al 
11995; .Berghofer et al.. .1996,) . .Moff^at & Michaud 



(119811) then further proposed that the inner wind re- 
gions are magnetically confined and forced to coro- 
tate with the star. Spectro-polarimetric observations 
did reveal an equato rial compression of the stel- 
lar wind of f Puppis dHarries & Howarth 199 6) but 
failed to indicate the pre sence of a strong magnetic 



field (ISchnerr et al.l l2008h . Shorter-term changes of 



the UV and optical line profiles with periods in the 
ran ge 2-20h were also reporte d by sev eral authors 
(B aad3 Il99ll iReid et all Il993l: iH owarth et all 1 1995 



IBerghofer et al]| 19961: iReid & Howarthlll996i and ref- 
erences therein) and usually attribut e d to n on-radial 
pulsations. Finally, lEversberg et al.l (Il998h detected 
stochastic substructures in the Hen/i4686 emission 
line moving away from the line center with time, 
which they interpreted as the signature of small-scale 
clumps. 

In the X-ray domain too, the variability of ^Puppis 



was in vestigated. Despite several early claims (Coll ura et al 
19891), Einstein observations did not reveal any sig- 
nificant variations: the false detections could be ex- 



1981 



pla ined by instrumental effects (iSnow et al 
IBerghofer etal.ilT996h. ROSAT da ta, however, told 
another story jBerghofer et al.ll 19961) . Variations of the 
mean count rate at high energies (0.9-2.0keV) were 
detected between the five ROSAT exposures, but only 
at a 2% significance level (i.e. lower than a 3cr detec- 
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tion). The longest ROSAT exposure (56.7ks spread 
over 11 days) was also analyzed to look for short- 
term variations. While and Kolmogorov-Smirnov 
methods did not yield definitive results, period folding 
techniques revealed the presence of a 1.44d"' fre- 
quency (=16. 7h period) in the ROSAT high-energy 
band (i.e. 0.9-2.0keV band). These high-energy vari- 
ations have a small amplitude (about 6%), hence the 
period is detected with only "95% confidence", but 
an additional fact reinforced the confidence in their 
presence: the X-ray varia tions appeared in phase with 
changes in the Ho' profile. iBerghofer et al. (1996) then 
proposed that the periodic modulation of the wind den- 
sity in the lower wind layers (as traced by Ha) triggers 
wind instabilities whic h, in turn, produce win d shocks 
and therefore X-rays. IBerghofer etalj(ll996l) however 
cautioned that their results should be checked with bet- 
ter quality X-ray data, and a later study of an AS CA 
dataset of 4"Puppis did not confirm the previous result 
( Oskinovaet al.l2001bh . Using the best X-ray dataset 



available at the present time (18 XMM-Newton expo- 
sures spread over a decade), we now try to investigate 
in detail the X-ray variability of ^Puppis. 

3. XMM-Newton observations and their reduction 

In the past decade of XMM-Newton observations, 
the star ^Puppis was observed during 18 exposures, 
mostly for calibration purposes. These datasets are 
excellent for studying the variability of ( Puppis since 
(1) the scheduled exposure times were often long (up 
to ~70ks) and (2) the observing dates probe weekly, 
monthly, and yearly timescales. Unfortunately, many 
observations were affected by soft proton background 
flares, resulting in total exposure times reduced by 
about 30%. A summary of the observations, as well 
as a detailed discussion of the reduction process (using 
SAS V 10.0.0) and of the spectra generation are given 
in Paper I. 

Up to five instruments aboard XMM-Newton simul- 
taneously record the X-ray emission of sources (three 
with low spectral resolution: EPIC-pn, EPIC-MOSl, 
EPIC-MOS2; and two with higher spectral resolution: 
RGSl, RGS2). Because of their separate optical paths 
and detector characteristics, they can be considered 
as five independent instruments observing simultane- 
ously. Indeed, the very different PSF observed with 
MOSl and MOS2 is a clear example that no two in- 
struments are identical aboard XMM-Newton, even if 
their names are similar. Furthermore, each observa- 



tion corresponds to a particular realization of the noise: 
therefore, it is possible (and even expected!) that each 
realization will slightly differ from another one. In the 
case of the XMM-Newton observations of ( Puppis, it 
means that a lightcurve recorded by one instrument 
will slightly differ from that recorded by another one, 
though they are compatible within the uncertainties. 
Having 5 independent datasets could be seen at first as 
a disadvantage, since any test or fitting will yield five 
not-exactly-identical answers, but that actually rein- 
forces the conclusions. Indeed, it reminds us that look- 
ing at a single observation with a single instrument, 
as so often is done except when working with XMM- 
Newton, does not allow any cross-checking and thus 
leaves room for spurious signal - in such cases, one 
must always be very careful not to over-interpret the 
data. Having multiple datasets enables us to check any 
detected signal, and in what follows, we emphasize 
only the secure variations, i.e., the significant changes 
observed by all instruments (EPIC+RGS). 

Our analysis is mainly based on the count rate 
lightcurves as they provide the most model-independent 
approach to the data. The EPIC lightcurves are com- 
posed of equivalent on-axis, full-PSF count rates, so 
that any offset or bad pixel problem was corrected 
to provide comparable lightcurves for all exposures 
(this possibility is not available for RGS). With this 
correction, all three EPIC detectors are considered to 
be stable over the spacecraft's lifetime by the XMM 
calibration team. Therefore, the flux found from spec- 
tral fitting and the count rate derived in the same en- 
ergy band are always related by a constant factor for a 
source with a constant spectral shape (which is the case 
of ^Puppis, see Paper I): using one or the other will 
yield exactly the same result, and fluxed lightcurves 
will thus not convey more information than those made 
from count rates. It should also be noted that fitting the 
spectrum of f Puppis is actually not an easy task (see 
Paper I), since we are reaching the limits of the instru- 
ments and atomic data precision. Within these limits 
and within the errors, the spectral shape of f Puppis 
is not changing from one observation to the next (see 
Paper I), but one could wonder whether this is true 
within an exposure. However, this would imply fit- 
ting spectra extracted in short time bins, hence very 
noisy. Such fitting would be highly unreliable, un- 
stable and it may even cause spurious changes in the 
spectral parameters, due to the complex spectral shape 
of f Puppis, the many local minima, and the uncertain- 
ties in atomic data. These spurious changes will thus 
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blur the actual source variations, rendering the task of 
studying them impossible in practice. Working with 
count rates avoid these problems, yielding clear and 
direct information on a given energy band, if not on a 
given spectral parameter other than flux. 

3.1. EPIC lightcurve production 

Lightcurves of the source and background were ex- 
tracted using the S AS task epiclccorr for four time bins 
(200s, 500s, Iks, 5ks) and in seven energy bands: total 
(0.3-4. keV), Berghofer's band (0.9-2. keV), soft (0.3- 
0.6 keV), medium (0.6-1.2 keV), hard (1.2-4. keV), 
very hard (4.-10. keV) and grand total (0.3-10. keV). 
The background and source regions are the same as for 
the spectral analysis (see Paper I). 

The choice of the energy bands results from a 
compromise, taking into account the appearance of 
f Puppis spectrum, count rates in each band, and pre- 
vious results. Indeed, one has to test bands for which 
results were reported before. As explained in Sect. 2, a 
periodicity was found using ROSAT in the 0.9-2 keV 
band. This is why we define a "Berghofer's band", to 
be able to confirm (or not) the presence of this signal in 
our much higher quality data. Another obvious choice 
is the total band, which enables us to use all available 
counts, and therefore get the smallest error, performing 
the most sensitive test of the overall level of the X-ray 
emission. Additional narrower bands (soft, medium, 
hard) are defined by taking into account the general 
spectral energy distribution (Fig. [U. The soft and 
medium bands probe the brighter part of the spectrum 
(with count rates above ~0.6 cts s ' keV '), while the 
hard band enables us to investigate the tail of the spec- 
trum. EPIC spectra showed a dim region around 0.6 
keV, hence the choice of this energy to split soft and 
medium energy bands. Count rates are rather similar 
for these "narrow" bands: 0.8, 0.4, 1, and 0.4 ctss"' 
for MOS (resp. 2, 2, 3.5, and 1 ctss ' for pn) in the 
Berghofer, soft, medium, and hard bands, respectively. 
The spectra appear very noisy above 4 keV, hence the 
choice of that energy as an upper limit for the main 
energy bands. Nevertheless, a (double) check was 
made as we want to ensure that including all data with 
>4 keV does not change the conclusions. Indeed, the 
grand total and total bands provide indistinguishable 
results. The very hard band would reveal the presence 
of very high energy phenomena (transient or not). For 
^Puppis, however, this band always presents a very 
low count rate, of the order of ~ (10 + 3) x 10"^ cts s"' 
for MOS (1 or 2), i.e. 3 orders of magnitude below the 



total band count rate, and similar to the background 
count rate in the same band for pn. This means that 
no bright emission at very high energies is present in 
this star, and we will therefore focus the analysis of 
^Puppis on the five main energy bands. 

The choice of the time bins was also made by tak- 
ing several factors into consideration. First, we consid- 
ered the problem of sampling different timescales. As 
one wished to study ^Puppis on as many timescales 
as possible, both short bins and longer bins are highly 
desirable. However, the upper limit on the time bin 
length is dictated by the length of the individual expo- 
sures (13ks for the shortest ones) since at least a few 
bins are required for a meaningful analysis within the 
exposure itself. Second, we looked at the impact of er- 
rors. Our goal is to detect the most subtle variations of 
f Puppis. To this aim, we need to combine the small- 
est possible errors. Considering the count rates for the 
diff'erent bands and instruments (see above), we can 
reach a 5% error on MOS data in 200s for the total 
band, in 500s for the medium and Berghofer bands, in 
Iks for the soft and hard bands; a 3% error on pn data 
in 200s for the total band, in 500s for the soft, medium 
and Berghofer bands, in Iks for the hard band. A time 
bin of 5ks has been added to get 1 % or less error on 
the total band and enhance the detection of longer-term 
changes within one exposure. 

To check the influence of the background (which is 
much fainter than the source in all observations), three 
sets of lightcurves were produced and analyzed indi- 
vidually: the raw source-i-background lightcurves, the 
background-subtracted lightcurves of the source and 
the lightcurves of the sole background region. The 
results found for the raw and background-subtracted 
lightcurves of the source are indistinguishable. An ex- 
ample of a lightcurve set is given in Fig. |2l 

Finally, two remarks must be made. To avoid 
very large errors and bad estimates of the count rates, 
we discarded bins displaying effective exposure time 
<50% of the time bin length. Our previous experience 
with XMM-Newton has shown us that including such 
bins degrades the results. It should also be noted that 
the SAS task epiclccorr provides equivalent on-axis, 
full PSF count rates, so that problems such as the pres- 
ence of an offset or of bad pixels were corrected to 
provide comparable lightcurves for all exposures. 
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3.2. RGS lightcurve production 

The task rgslccorr yielded raw source+ background 
lightcurves, background-corrected lightcurves for the 
source and background Ughtcurves for ten wave- 
length bins: total (6-30A, equivalent to 0.4-2 keV), 
Berghofer's band (6-14A, equivalent to 0.9-2 keV), 
soft (20-30A, equivalent to 0.4-0.6 keV), medium (14- 
20A, equivalent to 0.6-0.9 keV), Nvi (28-30A), Nvn 
(24.3-25.5A), Ovn (21.3-22A), Oviii (18.75-19.2A), 
Neix (13-14A), and Nex (11.8-12.5A). They were cal- 
culated for each instrument (with both orders com- 
bined) and for the whole RGS (with both instruments 
and both orders combined). Note that the wavelength 
of the Ovii line is at a gap of RGS2 while Nex is in a 
gap of RGS 1 (order 1). 

As for EPIC, there were obvious choices, such as 
the total band or bands covering the strongest lines, but 
also choices dictated by previous reports (Berghofer's 
band), and compromises (soft and medium bands - the 
soft band is similar to its "twin", the EPIC soft band, 
while the medium band appears in between soft and 
Berghofer's bands). The limits of the wavelength in- 
tervals were carefully chosen taking into account the 
global aspect of the high-resolution spectra (Fig. [3]). 
Limits of the bands fall in regions as free of lines as 
possible; the overall lower and upper limits are dic- 
tated by the increasing noise towards short and long 
wavelengths. It may be noted that count rates for the 
soft, medium and Berghofer bands are similar (about 
0.5 cts s"'). The bands covering lines were chosen to 
enclose the brightest isolated lines without enclosing 
too much of the neighbouring continuum. 

Again, as for EPIC, the choice of time bins repre- 
sents a compromise between getting a sufficient num- 
ber of counts to perform an analysis with a few per- 
cent relative error only and having enough bins within 
each exposure (of length 30 ks for the shortest ones). 
Bins of 500s and 2ks yield errors in the total band, for 
each RGS, of 7% and 3.5%, respectively; bins of 2ks 
and 5ks give errors in the soft, medium and Berghofer 
bands of 7% and 3.5%, respectively; bins of 5ks and 
lOks ensure errors of 10% and 5%, respectively, for 
bands linked to specific lines. 

Note that, contrary to epiclccorr, the SAS task 
rgslccorr does not correct for lower recorded fluxes if 
the source appears off-axis. While this does not appear 
to be a problem for small offsets (<1.5'), the count rate 
of ^ Puppis during Rev. 073 1 (where the offset reaches 
nearly 6') is clearly reduced, by about 25% in the to- 



tal band. In addition, the count rate from Rev. 0091 
at highest energies (shortest wavelengths) appears too 
high by about 20%. The origin of this problem is un- 
known (poor calibration at the earliest times of XMM- 
Newton operations or true brightening? - there is no 
EPIC data to confirm what happened). To avoid any 
interpretation problems, both datasets were discarded 
from global analyses (see Sect. 4 and 5). 

It is also important to remember that rgslccorr per- 
forms a randomization of the time tags, i.e. all runs 
of the same task, with the same parameters and input 
files will produce slightly different results (the differ- 
ence remaining within the error, of course). This is not 
the case of epiclccorr since a randomization is applied 
during the initial processing of the raw files (tasks ep- 
proc, emproc). 

4. Long-term variations (months to years) 

4.1. Count rates 

To search for the long-term variations on timescales 
of months to years, we computed the average count 
rates for each exposure and checked their constancy 
using tests. We remind the reader of a few cautions. 
For RGS, the data points associated with Rev. 073 1 are 
discarded from global analyses since the effective area 
changes due to the off-axis position of ^ Puppis are not 
taken into account in rgslccorr. The data from Rev. 
0091 are also discarded since they yield suspiciously 
high count rates (see XMM user handboolfl and Sect. 
3). For EPIC, the pn data taken with the Medium filter 
are most probably still slightly affected by pile-up (see 
also Paper I). Indeed, they yield clearly erratic results: 
while Revs. 535, 538 and 542 provide very similar 
count rates for their Thick filter data, this is not the 
case for their Medium filter data. 

Thereby considering only the best (i.e. most reli- 
able) data, it clearly appears that the count rates are 
not stable over timescales of years (Figs. |4] and |5]). 
The lightcurves obtained by each detector decline. The 
rate of decline however differs among them. Over 
the ~3800 days covered by the datasets, the EPIC- 
pn count rate decreased by about 6% in the total, soft 
and medium energy bands, 2% for the hard band and 
4% for Berghofer's band. The count rates of EPIC- 
MOS decreased by about 10% in the total, medium 
and Berghofer's bands, 6% in the hard band and 15% 
in the soft band. The RGS count rates decreased by 



http://xmm.esac.esa.int/extemal/xmm_user_support/documentation/uhb/index.html 
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18% in the total band, 28% in the soft band, 16% in 
the medium band, and 12% inBerghofer'sband. These 
variations are significant at the < 1 % levei as the error 
bars on the average count rates are very small. 

The fitting of EPIC spectra also showed a small 
decrease in flux of a few percent (Paper 1). More- 
over, the supernova remnant lEO 102-72, observed in- 
dependently from ^Puppis (though also for calibra- 
tion purposes), shows a decrease in its MOS count 
rate and flux similar to what is seen for ^ Puppis (M. 
Guainazzi, private communication). The exact origin 
of this systematic trend is not known (investigations 
are under way by the calibration teams), but such long- 
term effects are reminiscent of detector sensitivity age- 
ing problems. Up to now, the pn and MOS are consid- 
ered to be stable by the XMM-Newton calibration team 
(M. Guainazzi, private communicationjl, but there are 
apparently remaining imperfections in the long-term 
calibration. Since the decline of the registered count 
rates is currently not taken into account in the XMM 
calibration, we removed this trend by hand for further 
analyses and conclude that it is instrumental in nature. 

Note that true variations of the source may be su- 
perimposed to these global decreasing trends. For ex- 
ample, ^Puppis appears brighter during Rev. 1620 for 
both count rates and fluxes (derived from spectral fits, 
see Paper 1). This is linked to the changes observed 
on intermediate timescales (see Sect. 5) - with trends 
on the course of days, the source may indeed appear 
brighter or fainter sometimes, depending on the obser- 
vational sampling of these changes. 

Finally, we also calculated the ratios between the 
RGS count rates associated to lines of the same ele- 
ments - Ovii/Oviii, Neix/Nex and Nvi/Nvii. Com- 
paring globally the 16 secure RGS exposures (see 
above and Sect. 3), there is no obvious correlation 
between the different elements: the Neix/Nex and 
Nvi/Nvii ratios do significantly vary, but only the 
N vi/Nvii is clearly decreasing. The N lines constitute 
the most widely separated pair: as the long-term trend 
is instrumental in origin, such a calibration problem 
could affect some wavelengths more than others and a 



^ A significance level of 1 % implies that there is only 1 % chance to 
get the observed deviation from the hypothesis by coincidence. 

^See also Sect 1.3.2 and Fig. 1.14 of the 

EPIC calibration status which can be found on 
|http://xmm.vilspa.esa.e s/docs/documents/CAL-TN-001 8.pdf as 
well as Sect. 5 and Fig. 5 of the XMM cross-calibration status 
available on http://xmm.vilspa.esa.es/docs/documents/CAL-TN- 
0052.ps.gz 



larger separation would then lead to larger differences 
in lines. The average values of these ratios, based only 
on count rates (i.e. not on dereddened fluxes) and ex- 
cluding Revs. 0091 and 0731, are 1.44+0.01 for ni- 
trogen (both RGS, both orders), 1.28+0.01 for oxygen 
(RGSl, both orders), and 1.86+0.01 for neon (RGS2, 
both orders) - the choice of instruments takes into ac- 
count the fact that RGS2 could not record any informa- 
tion on Ovii, whereas the first order of RGS 1 contains 
no data on Nex (see Sect. 3). 

4.2. X-ray lines 

High-resolution X-ray spectra of good quality can 
be derived for each exposure, and used to search for 
variations. The lines recorded for each observations 
can first be compared to the average spectrum obtained 
when combining all observations (Fig. |6]l. No large, 
significant variation is detected. Small changes in flux 
(of the order of one sigma) can be spotted from time 
to time. The y are similar to those recorded at optica l 
wavelengths (lEversberg et al.ll998l;[Rauw et al.l2010l) . 
with double peaks where the blue to red ratio slightly 
changes. Their small amplitude and the low signal-to- 
noise of the data however prevent us from a detailed 
analysis: this must await the advent of more sensitive 
observatories. 

To be more quantitative, a temporal variance spec- 
trum analysis (TVS, Fullerton et al. 1996) was per- 
formed on the combined and fluxed RGS spectra (out- 
put of rgsfluxer). The TVS computes the squared 
difference between individual spectra and the aver- 
age one, taking the signal-to-noise ratios into account, 
which allows to detect statistically significant devia- 
tions from the average. Note that the spectra are cal- 
ibrated both in flux and wavelength within rgsfluxer 
(see Paper 1 for details), so that the pointing problems 
have no impact here. The relative weights for the dif- 
ferent spectra in the TVS were chosen equal to the 
count rate errors on the total band. The TVS appears 
overall flat (implying constancy), with only few dis- 
tinctive features (Fig. |7]i. For example, some peaks are 
found in the 20-24A region. This region is affected 
by gaps and missing CCDs, and there are therefore 
undefined values for wavelengths where no exposure 
exists. The exact position (in wavelength) where that 
occurs varies from one exposure to the next, notably 
because of pointing differences, and that results in ap- 
parent variability in the spectral set (hence peaks in the 
TVS). Increased noise at the shortest wavelengths (be- 
low 7A) and longest wavelengths (above 24A) is also 
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producing a larger TVS, unrelated to intrinsic variabil- 
ity of the source (though the peak at 7. 7 A remains un- 
explained). Globally, the TVS thus indicates that the 
lines in the X-ray spectrum of ^ Puppis are not signif- 
icantly varying from one exposure to the next. Very 
small changes are not excluded, however. Indeed, two 
small peaks of the TVS occur at 13.5 and 15 A, i.e. 
at the wavelength of the strongest lines: this indicates 
that the line profile variations of ^Puppis are just be- 
yond the reach of current facilities, and they may thus 
be detected in the future with the better sensitivity of 
new X-ray observatories. 

5. Intermediate-term variations (day to months) 

In this section, we investigate the data for the pres- 
ence of variations with timescales of day to months. 
Within each exposure, the same set of tests was ap- 
plied to each lightcurve. We first performed a x' test 
on all available individual lightcurves (i.e., 4 time bins 
and 7 energy bands for EPIC, 10 energy bins and 2 
time bins for RGS, see Sect. 3) for several null hy- 
potheses: constancy, linear variation, quadratic vari- 
ation. We further compared the improvement of the 

when increasing the number of parameters in the 
model (e.g. linear tr end vs constanc y) thanks to F-tests 
(see Sect. 12.2.5 in Lindgrenlll968h. A varia bilitv test 
using Bayesian blocks (BBs, Sc argle was also 

performed, through the FTOOLS battblocks. It was 
made on the 200s full (i.e., including bins with <50% 
eff'ective exposure time) EPIC lightcurves in the total 
energy band - those are the data with the largest signal- 
to-noise. For the most varying cases, a check was 
made by testing the event arrival times of the source, 
though this does not take into account the background, 
non-uniformities, or bad time intervals: results were 
similar to those derived from binned lightcurves. Un- 
less otherwise stated the adopted critical significance 
level is 1% throughout this paper. The individual re- 
sults for each exposure are detailed in the Appendix, 
and we only summarize them below. 

In general, the background is found to be variable. 
This is expected since many observations were af- 
fected by flares and some variations remain, despite 
the fact that the largest and narrowest ones have been 
cut out during the processing (see Paper I and Ap- 
pendix for details). 

For the source itself, however, only one thing is ob- 
vious: ^ Puppis does not vary much. For example, BBs 
are particularly useful to detect bursts, but there are 



none in ^ Puppis, and a single block is found in most 
cases to be the best representation of the lightcurves. 
The absence of bursts is confirmed by the;^'^ analyses. 

The error bars on each time bin are smaller for 
longer bins: small variations will thus be most eas- 
ily detected for long bins. On the other hand, varia- 
tions with timescales much smaller than the exposure 
duration will be smoothed out when using long time 
bins. Therefore, if short-term variations were dom- 
inating the overall variability of ^Puppis, the lowest 
variability level should be found for long time bins. 
What we find is the opposite trend: for both EPIC and 
RGS data, the longest time bins often yield the most 
variable lightcurves (see Fig. [8]l - though they rarely 
reach a significance level of 1%. This indicates that 
trends with timescales similar or larger than the expo- 
sure duration are more common in f Puppis than very 
short-term events. 

In fact, only 6 exposures show a significant non- 
constancy of their (mostly EPIC) lightcurves: Revs 
0795, 1164, 1343, 1620, 1814, and 1983 (Fig. |8]i. 
Moreover, in eight cases (Revs 0156, 1071, and the six 
previously quoted), modelling by a trend (linear and/or 
quadratic) yields a significant improvement of the 
over constancy, even if the significance levels attached 
to the individual x^ (for the constancy hypothesis) are 
not < 1 % in the additional two cases. It should be noted 
that 5 out of these datasets are amongst the longest ex- 
posures (see Paper I). Indeed, if we cut the data of Rev. 
1983 to keep only the first 10, 20, 30, 40, 50 or 60ks, 
the x^ (for the constancy hypothesis) never reaches a 
significance level <5% and there is no significant im- 
provement by fitting a line of non-zero slope rather 
than a constant! It thus seems that ^Puppis appears 
variable each time one looks at it for a sufficiently long 
time. For the most varying cases, i.e. Revs. 1343, 
1620, 1814, and 1983, a splitting of the Hghtcurves 
into two or three BBs is found (see e.g. Fig. |9]l. Again, 
this favors the existence of longer-term shallow trends 
over that of shorter-term "bursting" events. 

The performed tests can also provide information 
on the energy at which variations occur, when they 
are detected. Usually, it is expected that the hard 
band would be the most often variable in O-stars, as 
highly variable phenomena such as colliding winds 
or magnetically-confined winds mostly produce hard 
X-rays. However, this seems not to be the case in 
^Puppis, as there is no coherent behaviour for this 
band - in some exposures, the hard band appears as the 
most variable (i.e. that with the largest dispersion or 
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X^), in others as the least variable. In contrast, the soft 
band has a more consistent behaviour: it is only de- 
tected as variable when all other bands are varying too; 
it can therefore be classified as the least often variable 
one. Within the observational limitations, the varia- 
tions in ^ Puppis therefore appear as globally affecting 
its spectrum, rather than being particular to a specific 
(narrow) energy band. 

5.1. Fourier analysis 

The characteristic timescale of these intermediate- 
term variations is difficult to assess on individual ex- 
posures using x~ tests or BBs. These first analyses can 
only conclude that the timescale is longer than the typ- 
ical exposure length (i.e., a day or more). Looking fur- 
ther at the individual lightcurves, no obvious oscilla- 
tion is detected by eye, despite the full coverage of the 
putative ROSAT 17h-period by individual exposures 
and the large improvement in quality over the ROSAT 
data. It thus seems that a sinusoidal variation of sev- 
eral hours and an ampli tude of a few percents, such as 
that detected bv Berghofer et alj ( ll996i) . is transient, at 
best. 

To give a more quantitative assessment of the vari- 
ability timescale, we have performed period searches 
on global lightcurves (created by putting each individual- 
exposure lightcurve after one another, keeping their 
individual time tags), corrected for the long-term in- 
strumental tr end described in S ect. 4. 1 . We used the 
algorithms of | Heck et al. |(ll985l see also correction by 



Gosset et al.ll200Ih . As in previous sections, we con- 
sidered only the best datasets (see Paper I), i.e. the 
15 MOS exposures and the 10 pn exposures taken in 
small window + Thick filter mode, as well as the 16 
RGS exposures (i.e. excluding Revs. 0091 and 0731). 
The small time bins were favored (200s for EPIC, 500s 
for the total RGS band, and 2ks for the other RGS 
bands) since the Fourier-period-search algorithms are 
more affected by a small number of points than by 
large individual errors. Note that a more general pe- 
riod search technique which attempts to fit a period 
simultaneously with its harmonics does not yield ad- 
ditional insights. 

The resulting periodograms are shown in Fig. [TO] 
along with the spectral window (i.e. the aliasing struc- 
ture appearing naturally because of the temporal sam- 
pling). The periodograms consist of narrow peaks 
(since the dataset covers more than 10 years) tightly 
disposed in much broader features (because there are 



only a few exposures, of max. 70ks, within this 10-year 
timescale), which hamper an accurate determination of 
any period. However, two features are apparent for all 
bands except the hard one: (1) the highest peak occurs 
at about 0.3-0.4 d"', (2) the right wing of this peak 
decreases much more slowly than that of the spectral 
window - some additional power is thus present around 
0.7-1. 3d-'. 

To assess the significance of these signals, we per- 
formed a formal iterative decomposition in frequen- 
cies - we here illustrate only the analysis for the pn 
total bandpass. A first extracted frequency at vi = 
0.365 d"' (period P = 2.1 Ad, semi-amplitude a = 
0.25 cts/s) is significant against the null-hypothesis of 
white noise as deduced from ad hoc simulations. The 
following frequencies are V2 ~ 0.89 d ' (P = 1.12d, 
a = 0.12 cts/s) and V3 ~ 0.4 d"' (P ~ 2.5d, a = 
0.08 cts/s ; the last values are rather uncertain because 
they depend on the selected binning). These peaks in 
the Fourier periodogram are still characterized by a 
significance level lower than 0.001. The white noise 
hypothesis is thus formally rejected. 

However, these frequencies do not correspond to 
the true content of the signal of the star For exam- 
ple, when phased with the 0.365 d ' frequency, indi- 
vidual lightcurves appear one after another, with no 
common phase interval. As explained above, some in- 
dividual runs exhibit shallow trends. The time spanned 
by these runs combined with their scarcity over the 
covered decade has an unfortunate consequence: the 
Fourier transform is always able to find a few frequen- 
cies that are combining the various trends in a very 
constructive way. However the plethoric presence of 
gaps in the time series and the concomitant large num- 
ber of degrees of freedom is responsible for the posi- 
tive combination more than the possible coherency of 
the signal. As a test, we detrended the individual runs 
with a linear function before merging them. We com- 
puted the periodogram for the new time series. All 
the previously reported candidate frequencies disap- 
pear from the list of frequencies. The dominant one, 
for the 200 s binning, is V4 = 1.42 d"' (P = 0.70d, 
a = 0.054 cts/s) which corresponds to a peak whose 
significance level is around 0.001. If we detrend the 
X-ray lightcurves over the individual runs with a sec- 
ond degree polynomial, we obtain a time series that ex- 
hibits no outstanding peak and is thus in perfect agree- 
ment with the white noise hypothesis. Therefore, al- 
though the white noise hypothesis is rejected, the al- 
ternative hypothesis to adopt remains unsure. The star 
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exhibits weak although significant variations at a time 
scale between 0.7 day and several days but the sam- 
pling does not allow us to conclude to the existence 
of a coherent, fully deterministic signal above Poisson 
noise. The signal could be dominantly stochastic with 
some coherency at short timescale (some kind of red 
nois^. It is also possible that these changes are tran- 
sient: they would then not appear as a strictly periodic 
signal in the global lightcurve. Another possibility is 
that any signal could be not strictly stationary, in phase 
or frequency, from one observing run to the next, as of- 
ten seen for chang es of the optical spectra in Oef stars 
(IRauw et al ]l20Q3l - it must be recalled that our target, 
^Puppis, belongs to this spectral type category!). 

In summary, it is certain that ^ Puppis displays vari- 
ations with relative amplitudes up to 10-15% on the 
timescales of days - not hours (i.e., there is no trace 
of the signal attributed to non-radial pulsations, see 
Reid & Howarth However, our dataset suf- 



fers from the misfit between the possible frequencies 
and the actual sampling. Therefore, despite its ex- 
ceptional quality, it is still insufficient to detect con- 
fidently evidence for the presence of a coherent sig- 
nal with daily timescales. A modulation of the X-ray 
flux with such daily timescale was reported for the 
O-type dwarf ^Oph, and potentially associated with 
the corotatin g interacting features fo und through UV 
observations dOskinova et aklEoOlbh . ^Puppis shows 
a 5 d modulation also attributed to corotating regions 
iMo ffat & Michaud 1981). New, continuous X-ray ob- 
servations of ^ Puppis would be required to better char- 
acterize the timescale of the detected variability, and 
to establish whether these variations can be associated 
with such features. 

5.2. Autocorrelation analysis 

To search for temporal links between the photon 
arrival times, the data corrected for the instrumental 
effect were also analyzed using autocorrelation meth- 
ods (E_deLsOTi & Krolik 1988). The resulting autocor- 
relation functions are shown in Fig. [TT] Note that 
the binned functions were normalized by the actual 
number of points used and, as usual, by the disper- 
sion around the mean (observed variance sj^). The for- 
mer normalization is needed as the individual exposure 
times range from 15 to 77ks (Paper I), so that more 
data points can be used for some time shifts. 



A slight positive correlation appears for the total 
and medium bands for time shifts below 20ks. A slight 
anticorrelation is then detected for shifts of 50-60ks, 
i.e. similar to the typical duration of an exposure. 
A similar behaviour, though with even smaller ampli- 
tude, is seen for the soft, hard and Berghofer energy 
bands. This slow evolution towards lower correlation 
values confirms the existence of weak coherent varia- 
tions with timescales of tens of ks. 

The correlation amplitudes are however small, 
hence not highly significant, statistically speaking. 
However, a strong positive correlation for small shifts 
is expected since the wind configuration remains the 
same during some time (about one flow time, which 
is 7?»/voo=5.8ks for ^Puppis, as yco=2 250km s~' an d 



j^a,=1 8.6Rf7), using the parameters of lOskinova et al 



20061 and references therein). For longer shifts, two 
behaviors are then possible. On the one hand, the time 
bins could be independent for shifts longer than a flow 
time, as in a fully stochastic wind. In this case, the cor- 
relation function would rapidly drop to zero. On the 
other hand, if the time bins are coherently linked by a 
shallow trend, the correlation function should slowly 
decrease. This decrease is indeed what is observed, but 
we would thus expect a correlation signal with larger 
amplitudes. 

In practice, however, even a perfectly correlated 
signal (i.e., with a correlation of 1) will be diluted by 
the observational Poisson noise. The expected correla- 
tion value will then be given by 



C = 1 



(1) 



While white noise has, on average, the same amplitude whatever the 
frequency, red noise displays stronger amplitudes at low frequencies 



where the observed variance sj^ comprises the noise 
and the coherent variations of the source, whereas 
only corresponds to the noise. For the pn lightcurve 
in the total band with 200s time bins (globally de- 
trended, i.e. corrected for the instrumental effect), 
is 0.089cts^s"^ and the variance due to Poisson 
noise is estimated to be around 0.049 cts^ s"^. The 
peak height C should then reach 0.45, which is in good 
agreement with the value observed in Fig. [TT]for small 
time shifts. This implies that the intrinsic correlation 
of the noiseless signal could be nearing one: the trends 
are thus real. 

5.3. X-ray lines 

As already mentioned in Sect. 3, RGS lightcurves 
were also extracted for the brightest isolated X-ray 



9 



lines, enabling to determine whether these lines vary 
over shorter timescales. The observed variations are 
small, generally within Icr error bars. The only sig- 
nificant and coherent results are found for Revs. 0156 
and 1814. For the former, a linear trend with a posi- 
tive slope provides a much better fit to the Ne ix line 
flux. For the latter, the Neix line appears variable 
while the O viii clearly increases. Note however that 
these changes are detected with more significance us- 
ing RGS2 data than with RGSl data. 

Since the individual count rates are available, we 
also calculated the ratios between lines of the same ele- 
ments - O vii/O VIII, Ne ix/Ne x and N vi/N vn (Fig. (121) 



These ratios r eflect the t emperature ('B lumenthal et al. 



1972: Waldr on & Cassine lli 2007). and any change in 
the ratios would therefore be linked to temperature 
variations. Formally, opacity variations may also play 
a role in changing these ratios. However, the latter 
are not as sensitive to temperature as to absorption: a 
change in temperature by 15-20% results in a doubling 
of these ratios (note that the temperature dependence is 
similar for all ratios), whereas change of the absorbing 
column of a smooth wind by a factor of 2 (from 0. 1 
to 0.2 X 10^^ cm"2) yields variations of 7.5%, 25% and 
35% in the Ne, O and N ratios, respectively. As for 
count rates, the ratios were tested using tests. For 
the individual exposures, only two features appear sig- 
nificant. For Rev. 0156, the Ne ix/Ne x is better fit by a 
Une with negative slope. Similar but shallower possi- 
ble trends are detected by eye for the two other ratios, 
but they are not formally significant. For Rev. 1814, 
the Ne ix/Ne x appears much better fit by a concave/U- 
shaped parabola, but no similar signal is seen for the 
two other ratios. Note that these two revolutions are 
amongst the variable ones (see above), and that the 
reported changes may be linked to slight changes in 
spectra (see Sect. 4.2), for which higher-sensitivity in- 
struments are needed for a detailed characterization. 

6. Short-term variations (hours) 

Once the long-term instrumental decrease and the 
daily trends are removed, what signal is left? In prin- 
ciple, this is where the variability due to embedded 
wind-shocks should appear. In what follows, we will 
first present the new observational results, then the 
theoretical predictions of ^Puppis variability, and fi- 
nally compare the two to derive constraints on the wind 
structure. 



6.1. Observed Ughtcurves 

Relative dispersions were calculated for each of the 
observed lightcurves. These lightcurves have two spe- 
cific features which do not exist in the models (see 
next subsection): a temporal binning and a binning 
over a range of energies. In principle, this some- 
what smoothes out any variability if it is present with 
timescales shorter than the time-bin length, or if it 
changes strongly with energy (e.g. the lightcurves at 

0. 3 and 0.4 keV being uncorrected). However, the 
Poisson noise inevitably impacts on the data and there- 
fore prevents us to detect low-level variability (i.e. a 
few percent) with time bins smaller than 200s (see last 
fine of Table [Til - this limit would only be changed if 
a more sensitive instrument was used. In addition, the 
emitting parcels contain plasma emitting over a range 
of energies, even if isothermal, and this implies some 
correlation over different ranges of energies. A com- 
parison with synthetic curves can thus provide useful 
insights into the structure of the wind. 

As mentioned below, the steps in the synthetic 
lightcurves correspond to different wind configura- 
tions. The dispersion calculations for the data were 
thus performed in the following way: (1) the origi- 
nal lightcurves of each exposure were first detrended 
using the best-fit linear trend derived from calcu- 
lations (see above) as the simple model (presented in 
next section) does not predict nor model such features, 

(2) the count rate values of the 200s lightcurves were 
then sampled at 5 and lOks intervals (corresponding 
to about one and two wind flow times, respectively), 

1. e. only considering every 25th or 50th value, and 

(3) the relative dispersions of these new, reconstructed 
lightcurves were finally evaluated using 



el. disp = [ .^YjiCRi -a-bx ti)^/{N - 2) 



I mean 



where a and b are found from the best-fit trend, see 
Sect. 3.2.1, and mean = ZiCRi/crf)/ Zd/crj). This 
enables us to try to reproduce the 'independent wind 
configurations' of the successive 'time steps' in the 
synthetic curves. Dispersions were also calculated for 
the full (i.e. considering all time bins) 200s, 500s, Iks 
and 5ks lightcurves (after detrending). When one com- 
pares the dispersions of the full 200s lightcurves to 
these of the reconstructed lightcurves, they appear very 
similar The maximum differences amounts to +4% in 
the worst cases, or smaller (<1% in absolute value) 
when the observations are longer. This is unsurprising 
since dispersion estimates on, e.g., 2 or 3 bins being 
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less precise than on e.g. 50 bins: when one samples 
the lightcurves every 25 or 50 steps, a particularly dis- 
crepant realization of the noise has more impact when 
there are less bins. The sparse sampling therefore does 
not change much the results, so Table [T] reports only 
the dispersions of the full lightcurves. This Table gives 
in first column the revolution number, in columns 2-6 
the dispersions measured for the main 5 energy bands 
and the 200s binning, and in the 7th column the num- 
ber of data points used; Columns 8-13 report similar 
values, but for the 5ks binning. The last line provides 
for comparison the expected Poisson error, relative to 
the mean (i.e. 1/ ^jctrate x bin length). Note that only 
the EPIC data were used here, as the noise of the RGS 
lightcurves is even larger and would more easily mask 
subtle variations. 

Looking at Table [T] one thing is obvious: despite 
the fact that the exposure-long trends mentioned above 
are often not perfectly linear, the measured relative dis- 
persions around the best-fit lines are close to that ex- 
pected on the sole basis of the Poisson noise. This 
means that the true short-term variability of ( Puppis - 
and its potential wavelength dependence - remains hid- 
den in the noise, implying a very small amplitude for 
these 'intrinsic' short-term changes. 

These results should now be translated into con- 
straints on the wind structure. Considering X-rays to 
be emitted by many hot parcels, one may naively think 
that it is sufficient to notice that a 1 % dispersion "nat- 
urally" corresponds to lO"* "emitters", using simple 
Poisson statistics. However, reality is not as simple: 
emitters suffer from different amounts of absorption 
depending on their location and, worse, the absorption 
may be clumped too. In both cases, the simple reason- 
ing completely fails to apply, and a dedicated model is 
thus needed. 

6.2. A simple wind model 

Embedded wind-shocks are considered as responsi- 
ble for the X-ray emission of O-type stars. In these ex- 
panding and unstable winds, fast material encounters 
slow-moving material, giving rise to zones of dense 
gas, and the mutual collisions of such dense features 
then give rise to the X-rays. The radiation hydro- 
dynamic models (e.g. I peldmeier et al. 1997allbh pre- 
dict that the collision, and thus the heating, occurs 
quickly (tens of seconds), and the cooling of the re- 
sulting hot plasma is also rather rapid, as the cooling 
time remains lower than the wind flow time for dis- 



tances of several tens of stellar radii. A strong X-ray 
variability (2 orders of magnitude in amplitude) was 
predicted. However, early ROSAT and Einstein ob- 
servations failed to show such a strong variabihty of 
X-ray fluxes from massive stars jBerghofer & Schmitt 
1994bl and references therein). To reconcile the re- 



sults of the early observations with wind-shock the- 
ory, it was suggested that lateral break-up of the shells 
can lead to the presence of many parcels, resulting in 



low va riabihty jCassinelli et al. 1983 ; Feldmeier et al. 



1997al) : therefore 2-D or 3-D models were needed. 



First attempts for 2D w i nd m odels were made by 



Dessart & Owockil (120031 l2005h but only isothermal 



winds (i.e. without X-ray generation) have so far been 
considered. Another example is the recent work of 



Cassinelli et al.l (l2008h on bowshocks, but the model is 



not yet self-consistent in terms of production and evo- 
lution for an ensemble of clumps. Unfortunately, no 
new, multi-D hydrodynamical model has thus tackled 
the problem of X-ray generation in stellar winds, and 
we are thus left with alternative modelling paths. 

In this context, lOskinova et al. I (1200 lal) considered 
a spherically symmetric smooth cool wind permeated 
with discrete zones of hot X-ray emitting gas. It 
was found that X-ray variabihty depends on the fre- 
quency with which hot zones are generated, and on 
the cool wind opacity for the X-rays. It was shown, 
that in such smooth cool wind, the variability in soft 
band is expected to be smaller than in the hard band. 
Os kinova et al.. (2004) further developed a 2-D wind 
model where not only hot parcels are discrete, but the 
cool absorbing wind can be clumped too. We apply 
this model in the present paper to investigate the re- 
sulting X-ray variability. We only briefly recall here 
its basic features. 

Our model was designed to correctly perform the 
radiative transfer of X-rays in stellar winds, and it 
also reproduces the features derived from RHD sim- 
ulations. In this model, the X-ray emission originates 
from discrete zones of hot gas randomly distributed in 
an X-ray production zone extending from 1.2 to 100 
Rt . The choice of such a large zone may at first seem 
surprising, since the prediction of the ID hydrodynam- 
ical modelling place the X-ray emitting regions within 
a few stellar radii of the stellar surface. However, a de- 
tailed global analysis of the high-resolution spectrum 
of ^Puppis (Herve et al., submitted) shows that the X- 
ray emission zone must actually extend up to ~ 857?, 
to reproduce the observed spectrum, hence our choice 
of a large emission zone. Note that the lower bound- 
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ary of the region was chosen taking into account the 
analyses of individual fir triplets, which place onset 
radii in that range (e.g. Waldron & Cassinelli 2007, - 
note that these analyses also consider large emitting 
zones, with emission formally extending to infinity). 
All emitting parcels of gas contain the same amount 
of matter. The line emission is powered by colli- 
sional excitation and therefore scales with the density- 
squared. The density of the wind is derived based on 
the stellar mass-loss rate: from the mass conservation, 
it is p(r) = M /[Anr^vir)] where M is the mass-loss 
rate and the wind velocity v(r) is Vco(l - Ralrf [with 
and Rq chosen so that the photospheric velocity 
v(/?,) - O.Olvoo], as commonly used for a massive 
star wind. During motion, each hot fragment expands 
according to the continuity equation. Hence, the in- 
trinsic unattenuated X-ray luminosity of each hot frag- 
ment scales as l/r^v(r). The probability to find a hot 
fragment in the radius interval [r, r-Hdr] scales with 
l/v(r), i.e. emitters are concentrated at inner radii, 
where the wind is slow. The random radial location 
of fragments is d etermined by von Neumann's rejec- 
tion method (e.g.^ Press et al. I[l99 2h and their angular 
distribution is also random, with a uniform distribution 
over the sphere. Absorption and emission are decou- 
pled: there is no self-absorption for the emitting mate- 
rial and no re-emission of X-rays after absorption. The 
model allows for further sophistications, i.e. emissiv- 
ity can have a different scaling with density and den- 
sity can also be a parameter However, these param- 
eters affect the shape of X-ray emission line profiles, 
not relative variability, so in present study, we use the 
simplest emissivity. 

Once produced, the X-ray emission propagates 
through the cooler stellar wind which can absorb it. 
The velocity of this cool wind is assumed to follow 
the same velocity relation as for the hot wind compo- 
nent (the so-called "beta-law", see above). The cool 
stellar wind is assumed to be either a smooth cool 
wind or a set of cool clumps (randomly distributed 
over the 1.5 to 316 /?* range), both cases having the 
same overall mass-loss rate and optical-depth. For the 
fra gmented wind model, we use the "cones" model 
of lOskinova et al.l (120041) . with a lateral extent of 1 
degree for the spherical absorbers. A sketch of the 
model geometry is shown in Fig. [T3] Random radii 
are generated for each cone using the l/v(r) probabil- 
ity and von Neumann's rejection method, again. The 
total mass of a homogeneous wind enclosed between 
two subsequent radii is considered to be swept up in 



a dense fragment with the same optical depth as the 
homogeneous wind material, so that the fragment lo- 
cation is given by — dr' lv{r')]l[^^ dr' /r'^v(r')] 
where a and b are two subsequent random radii in the 
set determined by von Neuman's method. Note that, 
depending on the set of radii, cool clumps do not have 
all the same density, some being optically-thick while 
others are not blocking much light. 

The wind parameters (Voo, M, abundances) are de- 
rived fro m a non-LTE atmosph ere model specific to 
^Puppis ( Oskinova et al. 2006f) . matching the opti- 
cal/UV spectrum of the star. All these fix the clump lo- 
cation, size, density and optical depth, leaving as only 
free parameter the number of hot and cool clumps. It 
may be noted that this model reproduces well the ob- 
served X-ray line pro files, despite its simplifications 
dOskinova et al.ll2004 . 

Our analysis of observational data did not reveal 
strong spectral trends in the variability. Therefore, 
for simplicity and clarity we simulate here only the 
monochromatic X-ray flux at a few selected represen- 
tative wavelengths: 6A (~2keV, midpoint of the total 
EPIC band), 14A (-0.9 keV, midpoint of the medium 
EPIC band), and 19A (-0.65 keV, midpoint of the total 
RGS band). 

Some comments should be made about the syn- 
thetic lightcurves. First, the calculated flux is monochro 
matic and in arbitrary units. Only relative disper- 
sions, for example, can be calculated and compared 
to the data. Second, these variability lightcurves 
are not, strictly speaking, function of time. The 
ti me-dependent radiative hydrodynamic simulations 
of F eldmeier et al. 1(11997^ show that on the timescale 
longer than the flow time (5.8 ks f or ^^Puppis, as 
Vco=2250km s"'and /?,=18.6iJo. see lOskinova et al, 



2006! and references therein), the wind structure is 



renewed and is independent of the previous wind con- 
figuration. Our model reproduces this situation: each 
realization of our 2-D stochastic wind model is inde- 
pendent from the previous one and each model run rep- 
resents the wind configuration on a timescale shorter 
than the cooling time. In other words, we model the 
wind at some arbitrary moment of time, and compute 
the X-ray emergent flux at this moment; the next data 
point is calculated for a randomly different wind con- 
figuration for both absorbing and emitting parcels. We 
do not follow the wind expansion (this will be the sub- 
ject of a paper by Oskinova et al., in prep). While this 
approach does not allow us to model the detailed time 
evolution of the X-ray flux on timescales of hours, it 
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allows to model the relative amplitude of the X-ray 
variability for long stretches of randomly distributed 
observations, as appropriate for this XMM-Newton 
observing campaign. Examples of such iightcurves' 
are shown in the right panel of Fig. [13] while relative 
dispersions are presented in Table |2] 

Several conclusions can be drawn from Table|2] 

• As the number of emitting or absorbing parcels 
increases, the variability decreases, as could be 
expected. However, it is important to note that 
dispersions of one percent are only found when 
these numbers approach 10^. 

• A smooth wind is less variable than an otherwise 
equivalent clumped wind. 

• As wavelength increases (or energy decreases), 
the relative dispersion for a clumped wind re- 
mains stable or slightly increases whereas it de- 
cr eases for a smooth win d (as already reported 
in lOskinova et al.ll2001al) . In a clumped wind. 



the variability is thus less energy-dependent 
than in a smooth wind - in the extreme limit- 
ing case of a wind consisting of only opaque X- 
ray clumps, no energy depe ndence is expected 
dOskinova et alj|2004ll2006l) . 



• If absorbing clumps are distributed over a 
smaller region {Rmax of 100 rather than 3167?,), 
then the variability decreases. This is the effect 
of a smaller radial separation between clumps. 

Comparing these theoretical predictions with the 
observational results (see previous subsection), we 
found that the number of emitting and absorbing 
parcels is huge. Indeed, even in the most favourable 
case of smooth cool wind absorption, which is the 
least variable case, more than 100000 hot X-ray emit- 
ting zones must be present and contribute to the X-ray 
emission so that the relative flux variations remain be- 
low 1%. This number further increases when the cool 
wind fragmentation is also included in the model. 

It must be underlined that no previous study has put 
direct constraints on the number of clumps in O-stars. 
Some studies exist, however, fo r evolv ed massive stars. 
For example. iLepine &Mofi^atl(ll999h report that "be- 
tween 10"^ and 10"* clumps in the line emission region 
are needed to account for the line profile variability 
of the WR stars" that they analyzed. The X-ray data 
of ^Puppis suggest an even larger number, though 



this should apply to a larger zon e than formation re- 
gions of optical lines. In addition. iDavies et al.l (l2007h 
explained the polarization level of LBVs by either a 
few massive, optically-thick clumps (ejection rate of 
< 0.1 clump per flow time) or many small, optically- 
thin clumps (ejection rate of > 10^ clumps per flow 
time), with the latter option usually favored in litera- 
ture. Considering that the X-ray emission region and 
the cool wind absorption region cover several tens of 
stellar radii, our conclusion appears consistent with 
Davies' result, despite the different nature of the ob- 
jects under consideration. 

Thus, the sensitive XMM observations reveal that 
the stellar winds of O-stars are highly structured 
on small scales. It remains to be seen whether the 
theory of stellar wind instability can explain such 
a high degree of fragmentation. Some 2D mod- 
els of the line driven instability (LDI) in isother- 
mal stellar winds "show that radially compressed 
shells that develop initially from the LDI are system- 
atically broken up by Rayleigh-Taylor or thin-shell 
instabilit ies as these structures are accelerated out- 
ward" (Dess art & Owockil 2003), hence producing a 
lot of small-scale 2D structures, but the same au- 
thors later found "lateral c oherence of wind structures" 
dDessart & Owockill200 ?) and the question of the size 
and number of clumps therefore re mains unsettled 



The hy drodynamical simulations of iFeldmeier et al. 

^I1QQ'7^^' fn tViic A'^\7 t\\f^ cr\l*^ r\n*3C th^it t^^r-Vli^H thi^ r\rrvK. 



1997 a. to this day the sole ones that tackled the prob- 
lem of X-ray generation, though in ID) predicted that 
only a few strong shocks, simultaneously present in 
massive star winds, are responsible for most of the X- 
ray emission, whereas our results strongly challenge 
this. The late ral break-up of the X-ray emitting shells, 
advocated by lFeldmeier et aL ( 1997ah for lowering the 
X-ray variability, cannot provide an explanation (see 
Table |2|i. Moreover, if clumping is induced by sub- 
surface convection, hydrodynamical stellar evolution 
codes indicate that the total number of clumps in O- 
stars should typically amount to 6 x 10^ - 6 x 10"^ 
(|Cantiello et alJl2009l) . Our data indicate an order of 
magnitude larger values, therefore prompting further 
investigation on the nature of radiatively-driven stellar 
winds. 

7. Summary and Conclusion 

We have analyzed an exceptional set of X-ray ob- 
servations of ^Puppis: to date, there exists no more 



The wind temperature was set at the effective temperature of the star. 
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sensitive dataset nor a dataset with a better (time and 
spectral) coverage for a massive star. 

Over the decade of observations of ^Puppis, a de- 
creasing trend is clearly seen in the count rates. As a 
comparison with the fluxes determined from spectral 
fits (see Paper I) and with other X-ray sources shows, 
this is mostly due to instrumental/calibration problems 
(probably the ageing detector, whose sensitivity de- 
creases with time), and is not yet taken into account 
in the data calibration process. 

Comparing the X-ray lines appearing in the 18 
available high-resolution spectra yields again no sig- 
nificant, true long-term changes. Some shallow, Icr 
line profile variations are however reminiscent of those 
seen in optical - more sensitive observations should 
confirm this, and pinpoint the timescales on which that 
occurs. 

Once the instrumental efi'ect is taken out, we do not 
detect flare-like bursts of X-ray emission nor short- 
term variations (<ld, like e.g. the stellar pulsations 
detected in the optical domain) in the individual ex- 
posures. However, we detect statistically-significant 
variations of the count rates of f Puppis on timescales 
of > Id. Indeed, the Ughtcurves with the longest time 
bins appear more variable than those with shorter bins, 
and the longest datasets are systematically found to be 
variable. The detected changes appear as shallow in- 
creasing (Revs. 0156, 1071, 1620, 1814, and 1983) or 
decreasing (Rev. 1343) trends, or a mix of both (Rev. 
1164). No clear dependence with energy is found: in 
particular, the hard band is not the most often variable 
one. This suggests that the bulk of the X-ray emission 
is affected, not only a high-energy emission tail hnked 
to phenomena such as magnetic confinement. Further- 
more, no evidence for a coherent, systematic periodic- 
ity is found either. These slow modulations cannot be 
explained by the embedded wind shock scenario but 
are consistent with the presence of large-scale, slowly- 
moving structures in the wind, which may, for exam- 
ple, result from corotating interaction regions in the 
wind. Clearly, more data, specifically covering the full 
rotation period, are needed to settle the question of the 
origin of such features. 

Once instrumental and daily trends are taken out, 
^Puppis shows a surprisingly low level of variability. 
This places stringent constraints on the wind structure. 
A wind variability model tailored to f Puppis has been 
undertaken. It shows that only a very large number of 
emitting and absorbing clumps (> 10^) are able to re- 
produce the observed lightcurves of f Puppis. This is 



the first time that such a limit has been placed for an 
O-type star. The stringent hmit on number of clumps 
that we established questions some results from exist- 
ing models of stellar winds (link clumping-convection, 
fragmentation level) and X-ray production (number of 
hot gas zones), but these models were not calculated in 
3D. Future model developments should be done, tak- 
ing into account our high dumpiness result. 
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Fig. 1. — Position of the EPIC bands compared to the 
spectra from Rev. 1983 (top green points; pn, bottom 
black and red points: MOS). 
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Fig. 2. — EPIC data registered for Rev. 1071 by the 
EPIC-pn instrument, for different bands and time bins. 
The background-corrected lightcurves of the source 
appear at the top of each panel, while the background 
lightcurves appear at the bottom. The x-axis represents 
elapsed time (in seconds) since the beginning of the 
observation. 
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Fig. 3. — Position of the RGS bands compared to the 
combined high-resolution spectrum (bottom: full with 
error bars, top: zoom on the 10-30A region). 
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Fig. 4. — Evolution of the average count rate with time. For MOS (left panel), filled and open circles correspond to 
MOSl and 2, respectively. For pn (right panel), the 5 points at the beginning of the dataset appearing above others 
were taken with the Medium filter rather than the Thick one, and they are more affected by pile-up (see Paper 1 for 
details). In all panels, the best-fit linear trends (excluding the discrepant observations, see text) are drawn - for MOS, 
only that of MOSl is shown for clarity. 
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Fig. 5. — Same as Fig. |4lfor the RGS data. The point 
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with Rev. 0731,forwhich^Puppiswas6' off-axis (see 
text for details). 
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Fig. 6. — Examples of "large" variations in X-ray lines 
for a few observations (identified by their revolution 
number in each panel). The individual spectrum, com- 
bining both instrument and both orders for a given rev- 
olution, is plotted with its errors and compared to the 
combined, full RGS spectrum at the same wavelength 
(red line) in each top panel. Their difference (in the 
sense individual spectrum minus combined spectrum) 
is given in each bottom panel. 
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dotted hne corresponds to a significance level of 1%. 
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Fig. 8. — A few examples of intermediate-term varia- 
tions in hghtcurves. The x-axis corresponds to elapsed 
time (in seconds) since the beginning of the observa- 
tion. 
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Fig. 9. — Evolution of the total count rate in Rev. 
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the comparison, the corresponding lightcurves with 
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elapsed time (in seconds) since the beginning of the 
MOSl observation. 
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Fig. 10. — Periodograms in different energy bands for the EPIC instruments and the combined RGS data. 
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Fig. 11. — Autocorrelation functions for the EPIC instruments and the combined RGS data. 
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RGS 1 and RGS2, divided by the full exposure time. It 
is thus a kind of average between the two instruments 
(i.e. [RGSl + RGS 2]/2At). However, when no signal 
is recorded in one instrument, like is the case for O vn 
in RGS2, the RGS data are not shown. 
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Fig. 13. — Top: A snapshot showing the wind structure 
in the model for a small number of parcels. Black dots 
represent the absorbing clumps (which, despite their 
color, are not totally optically-thick), while crosses 
correspond to the emitting regions. Note that the size 
increases with distance to the star, as the parcels keep 
their angular size, as seen from ^Puppis. Only a small 
number of clumps is shown, for clarity. Bottom: Ex- 
cerpts of synthetic variability curves at 14A for a set of 
wind configurations. Next to each synthetic lightcurve 
are indicated the numbers of emitting zones followed 
by the number of absorbing clumps. For clarity, the 
lightcurves have been vertically shifted. 
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Table 1 : Relative dispersions (in %) measured for the 
observed XMM-Newton lightcurves after detrending 
(see text for detail). For the pn, the first five lines cor- 
respond to the data taken with the Medium filter (i.e. 
influenced by pile-up). 
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Table 2: Relative dispersions measured for a set of 
synthetic lightcurves. 
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" In all other runs the smooth wind starts after 316/?,, 
while in this run the cool fragment zone of the wind 
ends at 100/?*, which decreases the variability. 
* In this run, the number of clumps in the radial 
direction is the same (400) as in the previous case but 
the lateral size of each clump is 0.5°. No strong 
impact on the variability is found. 
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A. Results for individual exposures 

In this appendix, we present the variability properties of each exposure. We also remind the reader about back- 
ground flares afl'ecting the exposure (though they have been cut out during the processing, see Paper I). 

A.l. Rev. 0091 

This dataset only comprises data from the RGS. There were narrow background flares scattered all over the exposure 
and a broad flare affected the mid-exposure. At a significance level of 1%, no trend or deviation from a constant is 
detected. Count rates may be higher towards the end of the exposure, but only at the 10-20% level. A feeble oscillation 
can be detected by eye in the Berghofer band, with a recurrence time of about lOks, but it does not reach an amplitude 
larger than that of the 1-cr error bars. The count rate in the medium band and with the longest time bin appears as the 
most variable (i.e. they show ax^ associated with the lowest significance level, but this significance level is not <1%). 

A.2. Rev. 0156 

A narrow background flare aff'ected the pn data at mid-exposure. This dataset shows a trend towards increasing 
count rates. Though less obvious in MOS2 and RGSl, this trend is clearly detected through the significant improve- 
ment of the;^f^ when using a linear rather than a constant fit. This trend is detected for the total, soft, and medium bands 
(EPIC), and for the medium band (RGS), i.e. it concerns photons with energies below 1 keV. For the EPIC-MOSl 
lightcurve in the total band, the increase rate is 1 .7±0.7xl0~* cts s~^ (i.e. a 3% increase over the ~40ks exposure). For 
EPIC, count rates in the total band and/or calculated using the longest time bin appear as the most variable; for RGS, 
the medium band data and/or the hghtcurves with the longest time bins are the most variable. 

A.3. Rev. 0535 

No soft proton flare is detected for this exposure. Data are compatible with a constant count rate, without any trend 
detected or any obvious oscillation. The pn observation is cut in two parts for this revolution and the next two (see 
Paper I): here, the second pn dataset (taken with Medium filter) appears more variable than the first part (taken with 
Thick filter), in all bands and for the smallest time bin, but this difference is not formally significant. For EPIC, the 
hard band lightcurves (esp. those obtained with the smallest time bin) are the most variable; on the contrary, for the 
other bands, the hghtcurves calculated using the longest time bins are the most variable. RGS data do not yield any 
coherent result as to which time bin and energy band is the most variable. 

A.4. Rev. 0538 

A large flare occurred near the end of the observation. Data are compatible with a constant, without any trend 

detected or any obvious oscillation. For EPIC, the hard band lightcurves (esp. those obtained with the longest time 
bin) as well as lightcurves taken with the extreme time bins (200s and 5ks) are the most variable; for RGS, the shortest 
time bins generally yield the most variable hghtcurves. 

A.5. Rev. 0542 

A large flare occurred near the end of the observation, it did not affect the MOS data. Data are compatible with a 
constant count rate, without any trend detected or any obvious oscillation. The first pn dataset appears more variable 
than the second part, in all bands for the Iks bin, but it is not highly significant. Berghofer's band appears as the most 
variable for both EPIC and RGS; the longest time bins yield the most variable RGS lightcurves. 

A.6. Rev. 0552 

Only pn data are available for EPIC. A large flare occurred during the second half of the observation. The data from 
pn and RGSl indicate a slight increase (linear trend better than a constant at the <5% level) of the count rates in the 
second half of the observation, but this is not corroborated by the RGS2. While there is no time bin or energy band 
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favoring the variability properties of the RGS data, the medium band data and lightcurves calculated using the longest 
time bins (except for the hard band) yield the most variable lightcurves in the EPIC-pn. 

A.7. Rev. 0636 

A large flare occurred at the end of the MOS datasets, or in the middle of pn and RGS datasets. Data are compatible 
with a constant, without any trend detected or any obvious oscillation. There is no time bin or energy band favoring 
variability. 

A.8. Rev. 0731 

A large flare occurred during the last third of the observation. Only pn data are available for EPIC. Data are 
compatible with a constant count rate, without any trend detected or any obvious oscillation. For RGS, the medium 
band lightcurve calculated with the 2ks bin appears as the most variable; for EPIC, the hard band data are the most 
variable except for the longest time bins. 

A.9. Rev. 0795 

A large flare occurred during the second half of the observation. The count rate recorded in the hard band is not 
compatible with a constant - it shows a shallow decreasing trend - but only at significance levels of 1-10% (i.e. not 
formally significant). For EPIC, the hard band lightcurve appears as the most variable, but conclusions are unclear for 
the time bins: MOS data show more variability in the smallest time bins, while pn data favor variability in the longest 
ones. For RGS, the soft band appears as the most variable. 

A.IO. Rev. 0903 

A large flare occurred during the last third of the RGS observations. While the EPIC data (esp. pn) are much 
improved when fitted by an increasing trend, the RGS data are more compatible with a decreasing trend. Both types 
of instruments agree, however, that the longest time bins generally yield the most variable lightcurves. 

A.ll. Rev. 0980 

A large flare occurred near the end of the observation. Data are compatible with a constant, without any trend 
detected or any obvious oscillation. For EPIC, the largest variability (though not significant) is seen for the total band 
data while the medium data appear as the most variable in RGS. 

A.12. Rev. 1071 

Only pn data are available for EPIC. There is no locaUzed flare for this exposure but the background smoothly 
increases in the second half of the observation. An increasing trend is detected for the pn, in the total and soft bands, 
it is also detected in RGS (esp. RGSl) in the total band, but with a lower significance level. For RGS, the medium 
band lightcurve appears as the most variable for the longest time bins; for EPIC, the largest variations are detected for 
shortest time bins in the hard band and the longest time bins in the total band. 

A.13. Rev. 1096 

A few narrow background flares are scattered over the exposure, and a larger one occurs at the end of the observa- 
tion. Data are compatible with a constant rate, without any trend detected or any obvious oscillation. For RGS, there 
is no time bin nor energy band favoring variability; for EPIC, the most variable data are found for the longest time bin 
in the hard band. 
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A.14. Rev. 1164 



A few narrow background flares are spread over the exposure, and a larger one occurs at the beginning of the 
observation. Only MOS data are available for EPIC. A quadratic fit strongly improves the^^ in the total and medium 
bands of the MOS (which are not compatible with a constant count rate, at the 1-10% level), and in the total band 
of the RGS. For both types of instruments, the longest time bins provide the most variable lightcurves, but the most 
variable energy bands are the total and soft bands for MOS and RGS, respectively. 

A.15. Rev. 1343 

A large flare occurred during the first third of the observations. Data in the medium, total and Berghofer's EPIC 
bands are incompatible with a constant and significantly better fitted by a Unear decrease. The non-constancy in the 
Berghofer's band is also detected in the RGS data, though with a lower significance level. For EPIC, the medium band 
data and lightcurves calculated using 5ks bins appear as the most variable; for RGS, Berghofer's band appears the 
most variable, especially at small time bins. 

A.16. Rev. 1620 

A large flare occurred near the end of the observation. Data in all EPIC bands but the hard one are incompatible 
with a constant and significantly better fitted by a linear increase. The trend in the total band is also detected in the 
RGS data. For EPIC, the total and medium band data appear as the most variable, especially for the longest time bins; 
for RGS, the longest time bins yield the most variable lightcurves for the total and Berghofer's bands. 

A.17. Rev. 1814 

A small flare occurred at the beginning of the RGS observation. Data in all EPIC bands but the hard one are 
incompatible with a constant and are significantly better fitted by a large linear increase or, even better, a quadratic 
increase. The non-constancy and the improvement by an increasing trend are also detected in the RGS data for the 
total and medium bands. For EPIC, the hard data appear as the least variable, while the lightcurves calculated with the 
longest time bins are the most variable. For RGS, the medium band data are the most variable. 

A.18. Rev. 1983 

A large soft proton flare occurred near the end of the RGS observation. Data in the total and medium EPIC bands 
are significantly better fitted by a linear, or even quadratic, increase, while the hard and Berghofer's bands are only 
better fitted by a quadratic increase. The increasing trend is also detected in the RGS data for the total band. An 
oscillation is visible in Berghofer's band on both RGS and EPIC, with a recurrence time of about 50ks. In general, the 
longest time bins yield the most variable lightcurves. 
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